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Abstract 

A set of four scalar conditions involving normal components of 
the fields D and B and their normal derivatives at a planar surface 
is introduced, among which different pairs can be chosen to represent 
possible boundary conditions for the electromagnetic fields. Four such 
pairs turn out to yield meaningful boundary conditions and their re- 
sponses for an incident plane wave at a planar boundary are studied. 
The theory is subsequently generalized to more general boundary sur- 
faces defined by a coordinate function. It is found that two of the pairs 
correspond to the PEC and PMC conditions while the other two cor- 
respond to a mixture of PEC and PMC conditions for fields polarized 
TE or TM with respect to the coordinate defining the surface. 

1 Introduction 

When representing electromagnetic field problems as boundary-value prob- 
lems the boundary conditions are generally defined in terms of fields tan- 
gential to the boundary surface. A typical example is that of impedance 
boundary conditions which can be expressed in the form [1, 2] 

n X (E -f, ■ (n X H)) = 0, n ■ f , = f , ■ n = 0, (1) 
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for some surface impedance dyadic which may have infinite components. 
Here, n is the unit vector normal to the boundary surface. Special cases for 
the impedance boundary are the perfect electric counductor (PEC) boundary, 

f ^ = ^ n X E = 0, (2) 
and the perfect magnetic conductor (PMC) boundary [3, 4], 

Zs ^ oo, =^ n X E = 0, (3) 

also known as high-impedance surface [5, 6, 7]. A generalization of these is 
the perfect electromagnetic conductor (PEMC) boundary defined by 

% = -Ln X T, =^ n X (ME + H) = 0, (4) 

where M is the PEMC admittance [8]. In fact, for M = and 1/M = (4) 
yields the respective PMC and PEC boundaries. 

A different set of boundary conditions involving normal components of 
the vectors D and B at the boundary surface, 

n • D = 0, n • B = (5) 

has been recently introduced in [9, 10] in conjunction with electromagnetic 
cloaking and, independently, by these authors [11, 12, 13, 14]. A boundary 
defined by the conditions (5) was dubbed DB boundary for brevity. The cor- 
responding conditions for the E and H vectors depend on the medium in front 
of the boundary. In this study we assume a simple isotropic medium with 
permittivity e and permeability whence (5) is equivalent to the conditions 

n • E = 0, n • H = 0. (6) 

A more recent study of literature has revealed that the DB conditions ei- 
ther in the form of (5) or (6) have been considered much earlier. In fact, 
in 1959 V.H. Rumsey discussed the uniqueness of a problem involving the 
conditions (6) as well as their realization in terms of the interface of a uniax- 
ially anisotropic medium [15]. The uniqueness and existence problems were 
considered more exactly in subsequent papers [16, 17, 18, 19, 20]. 

The DB-boundary conditions (5) were introduced by these authors in 
[13] as following from the interface conditions for the half space z < of a 
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certain exotic material called uniaxial IB or skewon-axion medium [21, 22]. 
A more general proof is given in Appendix 1. A simpler realization for the 
DB boundary appears possible in terms of the planar interface 2; = of a 
uniaxially anisotropic medium defined by the permittivity and permeability 
dyadics [15, 13] 

I = e^u^u^ + eJt, ^ = /x^u^u^ + //Jt, (7) 
with the transverse unit dyadic defined by 

L = I - U^U^ = U^Ux + UyUy. (8) 

Assuming vanishing axial parameters, 0, /i^ —>■ 0, the fields in the 

uniaxial medium z < satisfy ^ and 0, whence because of the 

continuity, the DB conditions (5) are valid at the interface z = 0. Such a 
uniaxial medium was dubbed zero axial parameter (ZAP) medium [23]. Ob- 
viously, the same principle applies for curved boundaries as well, when the 
medium is locally uniaxial with vanishing normal components of permittivity 
and permeability. Zero-valued electromagnetic parameters and their appli- 
cations have been studied recently together with their realizations in terms 
of metamaterials [24, 25, 26, 27]. 

Basic properties of the DB boundary for electromagnetic fields have been 
recently studied. In [13] it was shown that, since the Poynting vector has 
only the normal component, the DB boundary is an isotropic soft surface in 
the definition of Kildal [28, 29]. Thus, coupling between aperture antennas 
on a DB plane is smaller than on a PEC plane. Further properties of a 
planar DB boundary were analyzed in [14]. It was shown that the DB plane 
can be replaced by a PEC plane for fields polarized TE^ with respect to the 
normal (z) direction and, by a PMC plane for TM^ fields. Thus, radiation 
from a current source J in front of the DB plane can be found by splitting 
the source in two parts, Jte radiating a TE^ field and Jtm radiating a TM^ 
field, whence the DB plane can be replaced by the images of the two source 
components, PEC image for the TE^ component and PMC image for the 
TM^ component. 

In [30] the resonator defined by a spherical DB boundary was studied. 
Splitting all modes in two sets, those polarized TE^ and TM'" with respect 
to the radial coordinate r, it was shown that the TE^ and the TM^ modes 
equal those of the respective PEC and PMC resonator. Thus, the number 
of modes with the same resonance frequency is double of that of the PEC 
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or PMC resonator which means that there is more freedom to define the 
resonance field. 

Finally, the circular waveguide defined by the DB boundary was analyzed 
in [31]. Since general modes cannot be decomposed in TE'' and TM'^ parts 
with respect to the polar radial coordinate p, the modes had to be computed 
in the classical way. It was shown that there may exist backward-wave modes 
in such a waveguide although there is no dispersion or periodic structure 
involved. 

The objective of the present paper is to extend the concept of DB bound- 
ary by adding another set of conditions involving normal derivatives of the 
normal field components. First, basic properties of the plane-wave refiected 
from planar boundaries satisfying different boundary conditions are derived, 
after which the more general curved boundary defined by a coordinate func- 
tion will be studied. 

2 Planar boundary conditions 

In the following we consider a planar boundary defined by 2; = and fields 
in the half space 2; > 0. The DB condition (5) can be expressed as vanishing 
of the z components of the two fields, 

= 0, B,^ 0, (9) 

at the boundary. If there are no sources at the boundary, the z components 
of the Maxwell equations yield conditions for the transverse components of 
the fields as 

Vt X Hi = 0, Vt X Et = 0. (10) 

Obviously, the conditions (10) are equivalent to (9). 

Let us now introduce another possible set of boundary conditions involv- 
ing the normal derivatives of the field components at the planar surface: 

d,D, = 0, d,B, = 0, (11) 

and let us call such a boundary by the name D'B' boundary for brevity. If 
there are no sources at the D'B' boundary, the fields D and B are solenoidal 
and they satisfy 



V • D = d,D, + Vt • Dt = 0, V • B = d,B, + V* • B^ = 0. (12) 
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Thus, the D'B' conditions (11) can be alternatively expressed in terms of the 
transverse components of the fields as 

Vf'Dt = 0, VfBt^ 0. (13) 

In an isotropic medium the D'B' conditions (11) and (13) are equivalent to 

d,E, - 0, d,H, = 0, (14) 

V* • Ei = 0, Vt • Hi = 0. (15) 

There are two other combinations of the four conditions in (9) and (11) 
which may appear useful. According to the previous pattern, let us call the 
conditions 

= 0, d,B, = (16) 
as those of the DB' boundary and the conditions 

d,D, = 0, = (17) 

as those of the D'B boundary. In an isotropic medium (16) can be replaced 

by 

Vt X Hi = 0, Vt • Ht = 0, (18) 
while the conditions (17) can be replaced by 

VfEt = 0, VtxEt = 0. (19) 

Prom (18) it follows that there exists a scalar potential ilj{x,y) in terms of 
which we can express the field at the DB' boundary and tjj satisfies the 
Laplace equation: 

Utix, y) = VMx, y), V?V(^, v) = 0. (20) 

Similarly, we can write for the field at the D'B boundary in terms of a 
potential as 

Et(x, y) = Vt0(x, y), V^0(x, y) = 0. (21) 

Assuming a localized source, the tangential components of the radiation fields 
are known to decay in the infinity as 1/r and the radial components a 1/r^ 
[35]. Following the argumentation of the Appendix, we must then have = 
and Ej = at the boundary surface. Thus, under the assumption of 
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localized sources the DB' conditions equal the PMC condition (3) and the 

D'B conditions equal the PEC condition (2). However, this is not valid for all 
non-localized sources. For example, a constant planar surface current gives 
rise to a TEM field with 0^ = and 3^ = everywhere, whence (16) and 
(17) are satisfied identically. This effect was discussed for the DB boundary 
in [14, 23]. 

3 Reflection of plane wave 
3.1 Plane- wave relations 

The basic problem associated with the D'B' boundary is to find the reflection 
of a plane wave from the planar boundary 2; = in the isotropic half space 
^ > defined by the parameters e, /i. Assuming exp(jcjt) time dependence 
and choosing the x axis so that the wave vector of the incident and reflected 
waves are in the xz plane, the field and wave vectors have the form 

E'(r) = E'e-^^'-", (r) = H'e'^^^-", (22) 

E''(r) = We-^^'^-", W{r) = We'^^'''', (23) 
k = — Uzkz, k = Ux^x + '^zkzi (24) 

Because the fields of a plane wave are divergenceless, they satisfy the orthog- 
onality relations 

k^ • E^ = 0, \i' W = 0, (25) 

k" ■ E" = 0, k" ■ = 0. (26) 

Assuming /c^ 7^ 0, i.e., excluding the normal incidence case, the fields can be 
expressed in terms of their z components as 

E^ = (u, + uMK)El + Uy(k/kx)vHi, (27) 

r/ff = (u, + UxkJkx)r]Hl - Uy{k/kx)El, (28) 
E'- = (u, - nxkz/kx)E: + Uy{k/kx)7]H:, (29) 
TjW = (u, - u,kz/kx)rjHl - Uy{k/kx)El. (30) 



Here, 



(31) 
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denote the respective wavenumber and wave-impedance quantities. 

Let us consider an incident plane wave whose z components satisfy the 
condition 

£;*sin^ + r//7*cos^ = 0. (32) 

Actually, any plane wave satisfies a condition of the form (32) for some 
parameter Q. (The TEM case E'^ = 0, i?^ = is excluded because of the 
assumption ^ 0.) The special cases of TE and TM waves correspond to 
the respective cases cos0 = and sin^ = 0. From (27), (28) the incident field 
vectors tangential to the boundary can be shown to satisfy the two conditions 

k 

EJ sin d + r]YL\ cos ^ = % — {j]E\ sin ^ - cos ^) , (33) 

kx 

Ej cos B - sin ^ = -u„ {r]R\ sin d - E\ cos Q) . (34) 

kx 

3.2 DB boundary 

The DB boundary conditions (9) at ^ = require 

El-rK^^. Hl + Hi^Q, (35) 

whence the refiected field components El, HI satisfy a relation of the same 
form (32) as the incident field components El, HI. 

From (29), (30) we see that the tangential components of the reflected 
field vectors satisfy 

E[ sin e + r/H[ cos ^ = -Uj, {riH'l sin 6 - El cos 6). (36) 

kx 

Combining this with (33) the condition for the total tangential field at the 
boundary becomes 

Etsin^ + r^HjCos^ = 0, (37) 

which is of the same form (32) as for the z components of the incident field. 
This leads us to the following conclusions: 

• For the TE incident wave with cos^ = 0, (37) yields E^ = which 
corresponds to the PEC condition (2). 

• For the TM incident wave with sin 6* = 0, (37) yields = which 
corresponds to the PMC condition (3). 
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• Denoting M = tanO/r], (37) yields ME^ + Hj = which corresponds 
to the PEMC condition (4). 

3.3 D'B' boundary 

At the D'B' boundary z — the fields satisfying the conditions (14) yield 

d^iEle^'^^' + Ele-^^^%=o = jUE^ - El) = 0, (38) 

d,{Hie^^^^ + H:e-^^^^),=o = jUK - HI) = 0. (39) 

whence, again, the reflected fleld components El, HI satisfy a relation of the 
same form (32) as the incident field components El, HI. 

Prom (29), (30) we now see that the tangential components of the refiected 
field vectors satisfy 

E[cos^-77H[sin^ = -u^^(77i/>in^-£;;cos^). (40) 

Combining with (34) the condition for the total tangential fields at the bound- 
ary becomes 

EtCOS^-r/Htsin^ = 0. (41) 
This condition leads us to the following conclusions: 

• For the TE incident wave with cos^ = 0, (41) yields = which 
corresponds to the PMC condition (3). 

• For the TM incident wave with sin^ = 0, (41) yields E^ = which 
corresponds to the PEC condition (2). 

• Denoting M' = — cot O/r], (41) yields M'Ef + Hj = which corresponds 
to the PEMC condition (4). 

Since these properties of the DB and D'B' boundary do not depend on 
the k vector of the incident plane- wave field (except that kx 0), they are 
valid for any plane waves. Being linear conditions, they are equally valid for 
general fields which do not have sources at the boundary. It is interesting to 
notice that the DB and D'B' conditions appear complementary in showing 
PEC and PMC properties to TE and TM polarized fields. Moreover, the two 
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PEMC admittances M and M' corresponding to the DB and D'B' boundary 
conditions satisfy the simple condition 

MM' = -l/r]'^. (42) 

This result has a close connection to the duality transformation [2] 

E ^ Erf = jr/H, H ^ Hrf = E/jr/, (43) 

which induces a transformation of media and boundary conditions. In par- 
ticular, the DB and D'B' conditions are invariant in the transformation but 
the PEMC admittance is transformed as M ^ = -l/rj^M = M' . 

3.4 DB' and D'B boundaries 

Let us finally consider the two other boundary conditions (16) and (17) for 
the incident plane wave (27), (28) satisfying the condition (32). For the DB' 
boundary the reflected fields satisfy 

^:+s; = o, hi-h: = q. (44) 

Prom (30) the reflected transverse magnetic field component becomes 

r;H[ = -u,{kJk,)rjHl + Uy{k/kM, (45) 

which compared with (28) can be seen to equal — ^^HJ. Because this is valid 
for any 6 DB' boundary condition equals the PMC condition (3) for any fields 
except when = 0. 

Similarly, the condition (17) leads to 

Ei-El^O, Hi + Hl^O, (46) 

and from (29) to 

El = -u,{kJk,)El - Uy{k/k,)rjHl, (47) 

which when compared with (27) equals — EJ. This corresponds to the PEC 
condition (2). 

As a summary wc can compare the four boundary conditions for TE and 
TM polarized incident fields in Table 1. 
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TE^ 


TM^ 


DB 


PEC 


PMC 


D'B' 


PMC 


PEC 


DB' 


PMC 


PMC 


D'B 


PEC 


PEC 



Table 1: Boundary conditions involving normal field components can be 
replaced by effective PEC and PMC conditions for fields with TE and TM 
polarizations. 

3.5 Reflection of TEM wave 

At this point one should consider the omitted case = corresponding 
to the TEM wave incident normally to the boundary. Since in this case 
the incident wave does not have normal field components, the preceding 
conditions arc already satisfied by the incident field and there appears to 
be no reflected field. However, since any reflected TEM wave can be added 
to the field, one may conclude that the conditions involving normal field 
components are not enough for defining the TEM fields. This anomaly, 
which was already pointed out by Rumsey [15], can be removed by adding 
the missing information. One may argue that a reflected TEM wave could 
also be associated to any TE and TM waves but, assuming a passive surface, 
it can be ruled out by power consideration. So the anomaly is associated 
with TEM waves, only, and we are free to define the reflection coefficient for 
the TEM wave (it may depend on the polarization of the wave). 

The additional information comes naturally if the surface is taken as the 
limiting case of a given physical interface. For example, for the DB bound- 
ary studied in terms of its realization by a ZAP-medium interface [14, 23, 32] 
the reflection of the TEM wave depends on the transverse medium param- 
eters et, /It of the medium. For the special case when the wave impedance 
of the TEM wave is equal in both media, there is no reflection. In this case 
we can deflne the corresponding DB boundary by requiring no reflection for 
the TEM wave. Another ZAP medium leads to another deflnition. Simi- 
lar consideration for the D'B '-boundary conditions requires a corresponding 
material reahzation which is yet to be found. For the DB' and D'B bound- 
aries it appears quite natural to assume that the respective PMC and PEC 
conditions are valid for the TEM wave as well. 

For a localized source giving rise to a continuous spectrum of plane waves. 
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the normally incident component has zero measure, i.e., it corresponds to zero 

portion of the total radiated power. Thus, the TEM component plays no role 
and can be neglected, as was already pointed out by Rumsey [15]. Also, in 
the case of a non-planar boundary an incident plane wave has a component 
appearing as a local TEM wave with zero energy which also can be omitted. 
It was shown through numerical computations for plane- wave scattering from 
spherical and cubic scatterers that there is no back-scattered wave when the 
DB or D'B' conditions are satisfied at the scatterer as predicted by the theory 
[33, 34]. 

3.6 Other possibilities 

For completeness, the two remaining possible combinations of boundary con- 
ditions, 

= 0, d,D, = 0, (48) 

= 0, d,B, = 0, (49) 

appear to be of little use. As an example, imposing (48) on the previous 
plane wave yields 

Ei + El^Q, Ei-El^O, ^ Ei^El^O. (50) 

This restricts the freedom of choice of the incident field. Thus, launching a 
TM incident field creates a contradiction, which can be so interpreted that 
(48) is of improper form. 

4 More general boundary surfaces 

Let us generalize the previous analysis by replacing the planar boundary 
surface by one defined by a function 2:3 (r) as 2:3 (r) = 0. Defining two other 
functions Xi{r) and 2:2 (r) so that they make a system of orthogonal coordi- 
nates satisfying 

Vx,(r) • Vx,(r) = 0, i^i, (51) 

allows us to express various differential operators in the form given in Ap- 
pendix 2. 

Expressing the fields as 

3 3 
E = J2^iEi, H = ^u,//,, (52) 
i=i i=i 
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and similarly for the D and B fields, the DB-boundary conditions are ex- 
pressed by 

Ds = 0, ^3 = 0, (53) 

for xs = 0. The conditions for the D'B' boundary are somewhat more com- 
plicated. Writing the expansion of the divergence 

V ■ F = Vi ■ + ^-^a,3(/ii/i2F3), (54) 

for the Cartesian coordinates with — z as 

V • F = Vt • Fi + d,F,, (55) 

gives us a reason to anticipate that a boundary condition of the form d^Fz — 
for the planar surface should take the form 9^3(/ii/i2-^3) = for the more 
general surface. 



4.1 Boundary conditions 

Assuming an isotropic medium bounded by the surface S : x^ — 0, we can 
propose four possible boundary conditions on S involving only the normal 
field components and and their normal derivatives as 

Ea^O and H^^O, (56) 

E3^0 and d^,{hih2H3) = 0, (57) 

d^^{hih2E3) = and H3 = 0, (58) 

9,3(/ii/i2£^3) = and 9,3(/ii/i2i^3) = 0. (59) 

Let us first consider consequences of these conditions. For £'3 = on the 
boundary S one of the Maxwell equations yields 

U3 • (V X H) = J-{d,,{h2H2) - d,,{hH,)) = 0, (60) 
H1H2 

while for H3 — the other Maxwell equation yields 

U3 • (V X E) = J-{d,^{h2E2) - d,,{h,E,)) = 0. (61) 

fllfl2 
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(60) is satisfied if there exists a scalar function ip{xi,X2) such that on S we 
can write ^ 

Hi{xi,X2) ^ —dxiip{xi,X2), (62) 
til 

H2{xi,X2) = ^d^2'^{xi,X2). (63) 

112 

Similarly, (61) is satisfied for a function (f){xi,X2) and 

Ei{xi,X2) = ^d^^(j){xi,X2), (64) 
til 

E2{xi,X2) = ^dx2(f){xi,X2), (65) 

tl2 

as can be easily checked. Thus, (60) and (61) can be compactly expressed in 
vector form as 

E^^O ^ Ilt{xi,X2)^Vtip{xi,X2), (66) 

Hs^O, ^ Et(xi,X2) = Vt(/>(xi,X2). (67) 

On the other hand, outside of sources, the divergence of the fields vanishes, 
whence from (111) we can write 

d,,{hih2Es{r)) = Q, ^ VfEi(r) = 0, (68) 

d,,{hih2H^{r))^0, VfHt(r) = 0. (69) 

Let us apply these on the different combinations of boundary conditions (56) 
- (59). 

4.2 DB boundary 

Starting from (56) (£'3 — O^H^ — 0) the fields on S can be expressed in 
terms of two potential functions as given in (66), (67). Let us now consider 
two special cases. Fields satisfying Es{r) = everywhere will be called TE^ 
fields and fields satisfying H^{r) = will be called TM^ fields. 

Obviously, a TE^ field satisfies 53(/ii/i2-E'3) =0 everywhere, including the 
boundary surface S. Prom (67) and (68) we conclude that on S the potential 
(f) satisfies 

Vt • Ei = Vt • (Vt0(xi, X2)) = Vl4){xi,X2) = 0. (70) 
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Prom the discussion in the Appendix we conclude that this imphes 

Et = VtH^,,X2) = 0, (71) 

i.e., PEC condition on the boundary surface -S". The result can be generalized 
to problems where S is not closed but extends to infinity, provided the sources 
are localized so that the fields vanish in the infinity. 

Similarly, for the TM'^ field the potential ijj must satisfy 

V,V(2:i,X2) = (72) 

on the surface S, whence the TM^ field sees the boundary defined by the DB 
conditions (56) as a PMC boundary. This condition can also be obtained 
from the duality transformation which swaps electric and magnetic fields 
and, hence, PEC and PMC conditions. The DB boundary is invariant to the 
duaity transformation. 

It must be pointed out that there is no guarantee that a given field can 
be expressed as a sum of partial fields TE and TM polarized with respect 
to the coordinate function xs{r) in the general case. Such a decomposition 
is known to be valid with respect to Euclidean coordinates and the radial 
spherical coordinate. 

4.3 D'B' boundary 

Considering now the boundary conditions (59), from (68), (69) the fields 
must satisfy 

Vt ■ Et{xi,X2) = 0, Vf iit{xi,X2) = (73) 

at the boundary. From (66) a TE^ field can be represented in terms of a 
potential iIj{xi,X2) which from (73) satisfies (72). From (115) wc again con- 
clude that a TM'^ field sees a D'B' boundary as a PMC boundary. Similarly a 
TM^ field sees the same boundary as a PEC boundary. Thus, in this respect 
DB and D'B' boundaries show complementary properties. 

It is clear that the D'B' conditions (59) represent a generalization of the 
conditions (14) for boundary surfaces more general than the planar surface. 
Comparing (59) and (111) whose last term can be written as V ■ (U3U3 ■ F), 
we see that the proper form for the D'B'-boundary conditions is 



V • (nn • E) = 0, V • (nn • H) = 0. 



(74) 
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As an example, for the spherical coordinates xi — 9, X2 — ip, X3 — r the 
metric coefficients are 

hi — he — r, h2 — h^ — r sin 9, h^ — hr — 1, (75) 

and the D'B'-boundary conditions (59) at the surface x^ = r = a have the 
form 

drir'^Er) = 0, drir^Hr) = (76) 
instead of drEj. — 0, d^Hr = as would be suggested by (14). 

4.4 DB' and D'B boundaries 

The mixed conditions (57) and (58) can be handled in the same way. In the 
case (57), (60) implies = VjV' on S while (69) implies V^V' = 0. Prom 
the reasoning given in the Appendix we obtain = on S. Thus, (57) 
corresponds to the PMC conditions for any fields. Similarly, we can show 
that (58) corresponds to the PEC conditions for any fields. 

4.5 The TF^ field 

The DB and D'B' conditions were tested above for the TE^ and TM"^ fields. 
Let us finally consider their generalization in terms of a combined field 

F = sin^E + cos^ryH, (77) 

where ^ is a parameter. Its component 3 is assumed to satisfy everywhere 
the condition 

F3 = sin OE^ + cos OrjH^ = 0. (78) 

Any field satisfying (78) is called a TF^ field. Since ds{hih2F^) — every- 
where, from (111) the TF^ field satisfies 

Vt • Ft = Vt • (sin 9Et + cos 9r]llt) = 0. (79) 

Inserting the Maxwell equations we can expand 

F = ^ V X (sin 9r}ll - cos ^E) , (80) 

whence the TF'^ field satisfies 

U3 • Vt X (sin 9rjHt - cos ^E*) = 0. (81) 
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Both (79) and (81) are valid in source-free regions. 

Because at the DB boundary Et and Hj are curl-free, (81) is automatically 
valid. The DB condition (56) implies 

U3 • Vt X (sin eEt + cos OrjUt) = (82) 

at the boundary. This combined with (79) and the reasoning for a closed 
surface given in the Appendix leads to 

sin eEt + cos erjUt = (83) 

at the DB boundary. This equals the PEMC condition 

MEt + r/Ht = 0, M = tan^/?]. (84) 

Following a similar path of reasoning, one can see that at the D'B' bound- 
ary (79) is automatically valid for a TF^ field while (81) and 

Vt • (sin eriUt - cos ^E^) = 0, (85) 

which follows from the D'B' conditions (59), correspond to the PEMC con- 
dition 

M'Et + rjUt = 0, M' = - cot (86) 

As a conclusion, for a TF^ field, which is a generalization of the TE^ and 
TM^ fields, both DB and D'B' boundaries appear as PEMC boundaries with 
the respective admittances M and M' . This includes as special cases the 
results for the TE^ and TM^ fields given above. 

5 Discussion 

The four boundary conditions (56) - (59) involving only field components 
normal to the boundary surface and their normal derivatives, form an inter- 
esting set. Two of these, (57) and (58) are alternatives for the respective 
PMC and PEC conditions which in terms of tangential fields are normally 
taken in the form (3) and (2). The other two conditions, dubbed as DB and 
D'B' conditions, (56), (59), appear new as kind of mixtures of PEC and PMC 
conditions or variants of the PEMC condition. The DB conditions have been 
introduced already in 1959 [15] but have not been applied for half a cen- 
tury; the D'B' conditions have not been discussed earlier to the knowledge 
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of these authors. Since their representations are so basic, DB and D'B' con- 
ditions seem to have their right of existence along with the PEC and PMC 
conditions. 

In introducing new boundary conditions the first question is to find their 
properties for the electromagnetic fields. This has been done here and in 
some of the previous publications. The next step would be to ponder about 
their possible applications. In case there are some of enough interest, the 
problem of practical realization of such boundary surfaces as an interface 
of some material comes up. As explained in the introduction, the planar 
DB boundary may have some application as an isotropic soft surface. Also, 
as shown in [9, 10], the DB conditions play a central role in the theory of 
electromagnetic cloaking. The DB boundary can be realized, e.g., by an 
interface of a uniaxial anisotropic medium with zero axial parameters (ZAP 
medium) or another medium called the uniaxial IB medium [13, 23]. Finding 
a realization for the D'B' medium is still an open problem. 

In [34] it has been shown that both DB and D'B' boundaries are self dual 
which implies that scatterers symmetric in 7r/2 rotation have no backscat- 
tering. This is an interesting property, shown to be valid by numerical com- 
putations, and it may have potential applications. As a continuation of the 
present study, it has been further shown that it is possible to define a gen- 
eralized DB boundary conditions of which the DB and D'B' conditions are 
special cases. Such boundaries are also self dual and have no backscattering 



In addition to the interesting physical properties of the novel boundary 
conditions, there may be some advantage in numerical electromagnetics by 
considering normal field components on the surface instead of the tangential 
components in the cases of DB' and D'B boundaries which correspond to the 
respective PMC and PEC conditions for the tangential fields. 

Appendix 1: Fields in uniaxial IB medium 

A medium called uniaxial IB medium has been defined in [13] by conditions 
of the form 



in terms of six parameters a, b, c, d, e, m. Let us briefly consider flelds in 
such a medium. Inserting (87) and (88) in the second one of the source-free 



[32]. 



D = aBt + hvL^B^ eu^ x E, 
H = mWz X B cEt -\- du^E^, 



(87) 
(88) 
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Maxwell equations 

V X E = -juB, (89) 

V X H = juB, (90) 

we can split both of them in two components. The axial components are, 
respectively, 

u, • Vt X = -jooB,, (91) 
mVt • Bt + cu, • Vt X Et = joubB,. (92) 
Eliminating the component yields the relation 

{b + c)u, • Vt X Et = -mVt • Bt (93) 

between the transverse fields. The transverse components of (89) and (90) 
can be written as 

d,u, xEt + jouBt = u, X WtE„ (94) 
{cd^ - juje)u^ X Ej - {md^ + jouajBt = du^ x Vt-B^. (95) 
Taking the divergence of (94) and (95) yields 

d,u, • Vt X Et = juVt ■ Bt, (96) 

{cdz - juje)u^ ■ X Ef = -{mdz + juja)Vt ■ Bt. (97) 

The three equations (93), (96) and (97) can now be reduced as follows. 
First we eliminate the term ■ Vt x Ej leaving us with the two equations 

[md, + Mb + c)] Vt • Bt = 0, (98) 

[mbd, + jLoime + a{b + c))]Vt ■ B< = 0. (99) 
As a second step we eliminate dtiVt ■ Bt) and arrive at the single equation 

[(6 + c) (6 - a) - me] Vt • Bj = 0. (100) 

Omitting the special case when the factor in brackets vanishes which corre- 
sponds to a certain subset of the medium, we conclude that the transverse 
B vector must be divergenceless, 

Vt • Bt = 0. (101) 
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Prom (93) it then follows that the E vector must satisfy 







(102) 



everywhere in the medium. Here we have also assumed that b+c ^ 0. Finally, 
from (91) and (87) we obtain conditions for the axial field components, 



To conclude, we have shown that any field in the uniaxial IB medium, 
defined by the medium equations (87), (88) in their general form, must satisfy 
the conditions (103). Thus, if we cut a layer, however thin, of the medium 
with two planes orthogonal to the z axis, the fields on both sides of the slab 
must satisfy the DB-boundary conditions (5) with n = u^. 

Appendix 2: Orthogonal coordinates 

The following review of differential operators for orthogonal coordinates de- 
fined by three functions a;i(r), a;2(r) and X3{r) has been taken from [35, 36]. 
The boundary surface S is defined by Xs{r) = 0. In this (r) and 

X2{r) define orthogonal coordinates on S. The coordinate unit vectors can 
be represented as 



where the hi — hi{r) are the metric-coefficient functions. These follow from 
the definition of the gradient. 



L>, = 0. 



(103) 



Ui = hiVxi{r), hi = (Vxi(r) • Vxi{r)) 



(104) 



3 



(105) 



The divergence and curl of a vector function 



3 



F(r) = ^u,F,(r) 



(106) 



1=1 



are defined by 




(107) 
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V X F = ^(d^^ihsF,) - d,,{h2F2)) 

+ ^{d.sihiF,)-d,,{hsF,)) 
mill 

+-^{d.,{h2F2) - a,,(/iiFi)), (108) 

fllfl2 

and the Laplacian by 



The differential operations can be spht in parts along the special coordi- 
nate Xs and transverse to it: 

V/^VJ + ^^a/, (110) 
'is 

V • F = Vt • + _L^a,3(/ii/i2F3), (111) 

111112113 

V X F = (V X F)t + ^{d,,{h2F2) - a,,(/iiFi)). (112) 

72 f V72 r , 1 o hih2 



If a scalar function (p{xi,X2) satisfies on the surface 5 : 2:3 = the Laplace 
equation 

V?0(xi,X2) = O, (114) 
and if the surface is closed, we can expand the surface integral as 

j \Vt(t>\^dS ^ J^f {(l>*Vt(l))dS - J (f)*V'^(f)dS. (115) 
s s s 

Both integrals on the right-hand side vanish. The last one because (114) and 
the middle one because the surface is closed. Thus, we can conclude that 
Vt0 = on 5". This elementary proof is, however, valid for simply connected 
surfaces only, and requires a more involved analysis for more complicated 
surfaces like the torus, see [19, 20]. 
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This condition remains valid for an open surface S extending to infinity 
when the integral over the boundary contour C of 5* in infinity, 

J Vt ■ {(l>*Vt(l>)dS = j^- {.(t>*yt(t>)dC, (116) 
s c 

vanishes. Here m is the unit vector normal to the contour and parallel to the 
surface in infinity. For a planar surface 2; = we have m = Up, the radial 
unit vector. As an example, if Vt0 is the tangential component of the electric 
field from a localized source, on the contour integral m ■ Vt0 corresponds to 
the radial component of the far field which is known to decay as 1/p^ along 
the plane. Since dC — pd(p, the integral vanishes and Vt0 = on the plane 
which corresponds to the PEC condition. Hovcver, this is not necessarily 
the case when the source extends to infinity. For example, for the normally 
incident TEM plane wave with constant Ej we have = Ej ■ p, whence the 
integral (116) becomes infinite. A more complete analysis of the open surface 
case remains still to be done. 
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